Unusual Microwave Response of Dirac Quasiparticles in Graphene 
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Recent experiments have proven that the quasiparticles in graphene obey a Dirac equation. Here 
we show that microwaves are an excellent probe of their unusual dynamics. When the chemical 
potential is small the intraband response can exhibit a cusp around zero frequency fl and this 
unusual lineshape changes to Drude-like by increasing the chemical potential with width also 
increasing linearly with ^. The interband contribution at T = is a constant independent of fl with a 
lower cutoff at 2/i. Distinctly different behavior occurs if interaction-induced phenomena in graphene 
cause an opening of a gap A. At large magnetic field B, the diagonal and Hall conductivities at 
small n become independent of B but remain nonzero and show structure associated with the lowest 
Landau level. This occurs because in the Dirac theory the energy of this level, _Eo = ±A, is field 
independent in sharp contrast to the conventional case. 
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The band structure of graphene (a single atomic layer 
of graphite) consists of two inequivalent pairs of the Dirac 
cones with apex at the hexagonal Brillouin zone cornersi. 
For zero chemical potential fj, one cone of the pair is full 
and the other empty. Through application of a gate 
voltage^ the upper unoccupied cone can be populated 
with electrons and by reversing the sign of the bias, holes 
can be introduced in the lower energy cone. Since a 
2+1 dimensional Dirac equation governs the dynam- 
ics of quasiparticles in graphene'^ many of its properties 
differ significantly from those of other materials. The 
most striking example is an unconventional quantization 
of the Hall conductivity predicted in Refsi^, unaware of 
the experiments reported later in Refs&i, which is re- 
lated to the anomalous properties of the lowest Landau 
level in Dirac theory. 

Up to date all experimental studies (see e.g. Refsi^) 
and most of the theoretical work on graphene have been 
made for dc transport properties. In this letter we study 
the microwave response of Dirac quasiparticles and show 
that it has several anomalous properties both in a strong 
magnetic field applied perpendicularly to the graphene 
plane and in the absence of the field. We argue that an 
experimental investigation of these features can provide 
added support to the existing consensus that the carri- 
ers in graphene are Dirac-like, but also shed light on not 
yet fully understood questions such as the character of 
impurity scattering and the presence of interaction in- 
duced phenomena. General expressions for the real part 
of the longitudinal (Jxxi^) and Hall axy{^) conductivity 
for graphene at any temperature T and magnetic field B 
are given in Ref>i£. They will not be repeated here al- 
though all numerical results presented are based on these 
exact expressions. To interpret these results we derive 
simple analytical expressions for axxi^) and <Txy{^) that 
allow us to elucidate the specifics of graphene. 

We begin with the simplest case of massless Dirac 
quasiparticles in zero field which corresponds to the con- 
tinuum limit of noninteracting quasiparticles on a hexag- 



onal lattice. The corresponding band structure is shown 
in Fig.n(a). 

In the limit of small impurity scattering rate r(w) 
neglecting the real part of the impurity self-energy, 
axxi^^T) takes a particularly simple form 
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where np{u}) = l/[e''^~'^^/"^ -f 1] is the Fermi distribu- 
tion (we set ks = ^) and Nj = 2 is the number of 
spin components. The first term of Eq. ^ describes 
the intraband transitions and the second one the inter- 
band. Similar expressions for the intraband contribu- 
tion have already appeared in the literature on d-wave 
superconductivit}*ii*i^ and for the interband case in the 
studies on d-density wavesi^ii^ii^. An essential feature of 
Eq. is that we kept the energy dependence of r(a;). In 
deriving this equation we have also assumed the small 
limit, so that T{ui') ~ r(u;) and for Vt <^T the difference 
[np{Lo) — npii^')]/^ can be replaced by the derivative 
~dnp{uj) / duo. 

Let us consider for a moment, only the intraband term 
of Eq. im 
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with (Too = e'^Nf /{2Tr'^h). If we take fi — and assume 
T(uj) — 7oo + a\uj\ with a small value of 700 as expected 
in Born approximation in the weak impurity scattering 
limit, we can get an approximate analytical expression 
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limit and T = when it reduces to (see Refii^) 
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FIG. 1: (Color online) Band structure of graphene. (a) The 
low-energy linear-dispersion -B(k) near the Dirac K and K' 
points for B = 0. (b) A possible modification of the quasi- 
particle spectrum by the finite gap (Dirac mass) A due to 
interaction-induced phenomena—*—. The chemical potential 
(indicated by horizontal line) jj, is shifted from zero by the 
gate voltage^, (c) Landau levels En in the Dirac theory of 
graphene. For a given direction of the magnetic field B ap- 
plied perpendicularly to graphene's plane the lowest (n = 0) 
Landau level has the energy i?o = — A at K and i5o = A at 
K' . The presence of a field independent n = level causes an 
unconventional Hall effect and anomalies in the strong field 
microwave response. 



Such an equation has been recently used in Ref»iS to ex- 
plain a cusp like behavior of the microwave conductivity 
observedi^ in high purity samples of ortho II YBCOe.s 
and some numerical results for graphene are presented in 
Ref.^^. This shows that if a cusp like behavior is observed 
in graphene, this would indicate that Born approxima- 
tion is relevant. If, however, the chemical potential is 
increased so that > T Eq. Q is replaced by 
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which has a Drude form. Here the factor of |^| gives 
the amount of spectral weight which is distributed ac- 
cording to a Drude form with scattering rate 700 + a|/^|. 
The width of the Drude peak can be increased simply by 



changing the gate voltage, fi oc sgnVg 



with impu- 



rity scattering left the same. What has happened is that 
T(lu) near uj = 11 rather than near uj — has now become 
the relevant quantity and it increases with \fi\. It is im- 
portant to note that the frequency dependence of r(ci;) in 
the more general case could be mapped out in this way, 
because the shape of r(a;) gets reflected in the shape of 
Cxx(S^)- For example, in the unitary limit of a T-matrix 
approximation the dependence of r(a;) is completely dif- 
ferent from Born limit and is inversely proportional to 
the electron density of states^. This dependence in turn 
modulates the dependence of the intraband term. 
Next we return to Eq. ^ and consider the strict F = 
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For fi — there is no intraband Drude contribution here 
proportional to S{n) because F = 0. Also the inter- 
band contribution extends from zero to the band edge 
(assumed very large in this work). It is independent of 
energy Q and in the bare bubble approximation used to 
derive Eqs. (QJ, 10) is equal to 7re^iV//(4/i) in heightiS,. 
The flatness of this response can be traced to the topol- 
ogy of the Dirac cone and the linear in momentum de- 
pendence of the Dirac quasiparticles. As fi increases, the 
interband transitions become gapped at = 2|/Lt| and the 
missing spectral weight reappears in the (5(ri)-function 
contribution. 

A nonlinear behavior of the diagonal magnetoresistiv- 
ity, pxx in Kish graphite in high magnetic field has been 
observed and interpreted in Ref^ in terms of the forma- 
tion of the new charge-density-wave phase. Also as pre- 
dicted in Refs^Siii, an interaction-induced phenomena 
may cause an opening of an excitonic gap (see Fig.^(b)) 
in the quasiparticle spectrum. Depending on the model 
parameters, the gap opens in zero or in a finite magnetic 
field. Recent measurements^i of the Hall resistivity made 
on highly oriented pyrolytic graphite indicate that there 
is a gap at i? = 0.03T. In our opinion, if the value of 
the gap is expected to be the order of a few Kelvin, mi- 
crowaves are the best technique to detect its existence. 
In this instance, for massive Dirac quasiparticles in zero 
field Eq. (jSJ acquires the form 



.(r^) = ^.(^7)(^^-^^y-^^) 



4A2 



4/1 



M n , AN 

— -max(|^|. A) 



(6) 



Here the thcta function 6* (|ri|/2 — max(|/i|. A)) cuts off 
the low r2 part of the interband contribution at 2|/i| 
or 2A whichever is the largest. The additional factor 
(ri^ -I- 4A^)/r2^ is to be noted and allows one to distin- 
guish in cr(ri) the case of finite A from finite /i. As we 
saw in relation to Eq. (O when A = the interband con- 
tribution remains flat although cuts off at 2|/u|. For finite 
A this is no longer that case and for < A the inter- 
band spectral weight at the gap has been increased by a 
factor of 2. This is a signature of a gap opening and ef- 
fectively represents a change in the geometry of the Dirac 
cones as illustrated in Fig. ^ (from (a) to (b) panel) . In 
Fig. 121 we show results for the microwave conductivity 
axx{^,T) vs VL at finite temperature. The long dashed 
curve for T = IK shows a Drude peak at small ft (intra) 
superimposed on a interband contribution of 7re^/(2/i) 
constant independent of O except for some depletion be- 
low 200GHz due to finite temperature with optical spec- 
tral weight transferred to the Drude peak. When T is 
increased to 5K (solid curve), the thermal depletion in- 
creased considerably with increased weight in the Drude 
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FIG. 2: (Color online) The microwave conductivity axxi^,T) 
in units e^/h vs microwave frequency il in GHz. The long 
dashed is for T — IK and the solid for T = 5K. In both cases 
fi — A = and T = O.IK. The dash-dotted curve has a gap 
A = 5K and r = IK. 



peak. On the other hand, if instead of increasing T we 
open a gap A = 5K (the dash dotted curve), the Drude 
peak is almost completely depleted as is the region below 
2A ~ 200GHz with extra optical weight accumulating at 
and above 2A in the interband region, where it is dis- 
tributed over the region of H > (2 4)A. The peak at 

= 2 A in axx{^) is characteristic of the opening of a 
gap as we have already described. 

Massless and massive Dirac quasiparticles in a strong 
field. The latest measurements in graphene samples-^ 
show that at B = 25T near /i = the diagonal and Hall 
dc magneto-resistances do not conform to the standard 
quantum Hall observation. Here we target the microwave 
response in the same region. When a magnetic field B is 
applied, the expressions for the conductivity become very 
complicated as given in Refii^. In the very high field limit 
the Landau levels with n ^ can all be shifted to high 
energy, but as we have already illustrated in Fig.n(c) this 
is not the case for n = level. This level may only depend 
on the field indirectly if the gap A{B,^) is induced by the 
field^''. (the phenomenon of magnetic catalysis^), so that 
it remains well below the rest of the levels. If we take 
B very large, only the contribution ^ 1/B associated 
with the transitions from n = to n = 1 levels survives 
and after multiplying it by the degeneracy factor of the 
Landau levels ^ B, we arrive at a very simple formula 
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where "0 is the digamma function. Here for simplicity 
we neglected the Zeeman splitting in the strong field. 
Eq. Q which is applicable for an arbitrary relationship 
between f2, T and F, shows that the diagonal conductiv- 



ity <7xx{^,T) is independent of B, but the position in 
energy of the two inequivalent n = Landau orbitals 
namely /i — A(i?, /i) and fi + A(i?, /i) measured relative 
to the chemical potential (see Fig. ^ (c)) remains in the 
problem. Results based on calculation of the full expres- 
sion for axx{^,T) from Refiiii are shown in Fig. Eland 
we have verified that they are in a good agreement with 
the approximate expression {T)). There are four curves. 
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FIG. 3: (Color online) Conductivity axx{^,T) in units e^/h 
vs microwave frequency Q in GHz. The chemical potential is 
set at 6K, the temperature at 0.5K and the impurity scatter- 
ing rate assumed constant set at 0.5K. Long dashed curve 
B = 0.05T and A = 0, sohd curve A = 0, dash-dotted 
A — 3K and short dashed A — 7K. All three last curves 
are for B = IT. 



Comparison of solid and long-dashed curves shows that 
the large field limit is already attained for B — 0.05T. 
Here the gap A = and the chemical potential equal to 
6K is clearly identifiable at 120GHz, where a large rise 
smeared by temperature (0.5 K) and impurities F = 0.5K 
is seen in axx- When a gap develops, however, the shape 
of the curve is completely changed. In both long-short 
dashed and short dashed curves, structure at — A and 
II + A is clearly visible reflecting that the energy of the 
n = Landau level at K point is A and at K' point 
is —A, respectively. This is yet another clear signature 
of the Dirac character of the quasiparticles in graphene 
which can be used to detect the gap. Indeed since in large 
fields all changes in the shape of (Jxx{^) are caused by 
the gap, one may consider that if the gap A{B) develops 
as the field increases, the curves with A = 0, A = 3K 
and A = 7K shown in Fig. O would actually show up as 
the field increases. 

Finally we consider the microwave Hall conductivity 
which can also provide more insight into the properties 
of graphene. It can be shown that under the same condi- 
tions as were used to derive Eq. Q, the microwave Hall 
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conductivity is given by 



''NfSgii{eB) 



Im 



(A 



-Re 



2n 

r + i{fi-n + A) 1 

2^ ^ 2 



T + i{fj. + A) 1 
27rT ^ 2 
r + i{n + n + A) ^ 1 
2 




In Fig.2]we show results for the frequency dependence of 
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FIG. 4: (Color online) Conductivity axy{0,,T) in units e^/h 
vs microwave frequency Q in GHz. The chemical potential 
is set at — lOK, the temperature at 0.5K and the impurity 
scattering rate assumed constant set at 4K. The gap A is 
also indicated in the figure. 

the Hall conductivity in units of e^/h. As found in Ref.^" 
the opening of a gap decreases the dc value of the Hall 



conductivity. For the case A = (dashed line) there 
is additionally structure at f2 ^ \fj,\, while for A finite 
the equivalent structure is moved to energies + A and 
— A (see Fig. ^ (c)) as we have seen in the case of 
diagonal conductivity. 

Up to now we have considered the influence of the gap 
opening on the conductivity when the gap itself was set 
to a constant. A more interesting possibility is to explore 
the dependence of the gap A on /x and B by changing the 
gate voltage or the magnetic field. It is expected that a 
decrease in ^ or an increase of B favors the opening of 
the gap when < ot B > Be as was predicted in 
Refjii. The important conclusion of our studies is that 
the generation of an excitonic gap should lead to a new 
insulator phase around the point = or in very large 
magnetic fields and this would be one more striking dif- 
ference of QHE in graphene from standard semiconduc- 
tors where the fractional quantum Hall effect is instead 
observed when entering the lowest Landau level. 

In summary the microwave response of graphene shows 
many distinct characteristics which reflect the Dirac na- 
ture of the quasiparticles. These include intra and in- 
terband contributions to the conductivity in zero mag- 
netic field and microwave response at large fields. While 
in conventional case large B would shift all Landau lev- 
els to high energies above the microwave region for the 
Dirac case the lowest Landau level contribution remains 
unshifted and affects a^xiS^) and axy{^) in the small f2 
region in a distinct way which can be used to study the 
effect of a gap opening induced by the magnetic field. 
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